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The neutrino role in primordial nucleosynthesis is reviewed. The importance 
of nonequilibrium effects is emphasized both for the standard massless and pos- 
sibly massive neutrinos. The upper bound on tau neutrino mass is presented. A 
spatial variation of primordial abundances and a possibility of observing them 
by precise measurements of the CMB anisotropy are considered. The nucleosyn- 
thesis bounds on the parameters of neutrino oscillations into sterile neutrinos 
are discussed. 
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^ ; 1 Introduction 

It is well known that neutrinos have a very important impact on cosmology (for a 
recent review and a list of references see e.g. In particular, the comparison 

of primordial nucleosynthesis theory with observational data permits to put rather 
stringent bounds on neutrino properties, their mass, number of flavors, possible new 
interactions, etc. To this end a detailed study of neutrino kinetics in the primeval 
plasma, especially nonequilibrium corrections, which happen to be quite essential, is 
of primary importance and the proper treatment of the latter significantly changes the 
results of the simpler equilibrium calculations. Technically the problem is very com- 
plicated and demands an accurate numerical solution of a system of coupled integro- 
differential kinetic equations, however in many cases a rough order-of-magnitude es- 
timate of non-equilibrium corrections can be done analytically. 
1 Also: ITEP, Bol. Chcrcmushkinskaya 25, Moscow 113259, Russia. 
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Surprisingly nonequilibrium corrections to the energy spectrum are not very small 
even for massless neutrinos in the standard model and of course they are very large in 
the case of possibly heavy v T with the mass in MeV range. Account of nonequilibrium 
effects permits to considerably improve the nucleosynthesis bounds on the mass of v T . 

Another important effect, where deviations from the standard equilibrium non- 
degenerate Fermi-Dirac distribution 

fM = eME/nt)] + l (1) 

is essential for nucleosynthesis, is a possible lepton asymmetry. The latter could be 
either primordial, generated at a very early stage, or it may arise during nucleosyn- 
thesis epoch due to nonequilibrium neutrino oscillations. Even if lepton asymmetry is 
not generated by oscillations, they still might have a very important impact on nucle- 
osynthesis and the study of nucleosynthesis leads to interesting bounds on oscillation 
parameters. 

If there are some new particles, abundant in the plasma during nucleosynthe- 
sis, they would change the cooling rate of the plasma and thus change the standard 
abundances of the light elements. Normally the effect of such new particles is just 
to change the expansion/cooling rates but in some cases their interaction may also 
produce nonequilibrium neutrinos, especially u e , and in this case the impact on nu- 
cleosynthesis would be significantly different. 

In what follows I briefly review these subjects. In section 2 nonequilibrium cor- 
rections to the spectra of normal massless neutrinos and their possible observational 
manifestations are discussed. In section 3 the role of a possibly heavy v T in nucleosyn- 
thesis is considered and an upper bound on its mass is presented. Neutrino degeneracy 
and especially a possible variation of the latter on the cosmological scales (a few 100 
Mpc or even Gpc) is discussed in section 4. In section 5 neutrino oscillations are 
considered. In Conclusion the main results of this brief review is summarized. 
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2 Nonequilibrium massless neutrinos 



It is usually assumed that thermal relics with m = are in perfect equilibrium state 
even after decoupling. For the photons in cosmic microwave background (CMB) it 
is known with a very high accuracy. The same assumption is made about neutrinos 
so that their distribution is given by eq. ([!]). Indeed when the interaction rate is 
high in comparison with the expansion rate, T int 3> H, the equilibrium is evidently 
established. When interactions can be neglected the distribution function may have 
an arbitrary form but for massless particles the equilibrium distribution is preserved 
if it was established earlier at a dense and hot stage when the interaction was fast. 
One can see that from kinetic equation in the expanding universe: 

(d t -Hpd p )f J (p„t)=lf l (2) 

where the collision integral in the r.h.s. vanishes for the equilibrium functions: 

f(e q) = ( e E/T-,/T ±i y 1 (3) 

The temperature T and chemical potential fi may be functions of time. 

The l.h.s. is annihilated by / = if the following condition is fulfilled for any 
value of particle energy E and momentum p = \/E 2 — m 2 : 

t + H *?Z-!L(*-t)=o (4) 

T E dp E \n T J v ; 

This can only be true if p = E (i.e m = 0), T/T = —H, and \l ~ T. It can be 
shown that for massless particles, which initially possessed equilibrium distribution, 
temperature and chemical potential indeed satisfy these requirements for I coU = 0, 
so the equilibrium distribution is not destroyed even when the interaction is switched 
off. 

It would be true for neutrinos if they instantly decoupled from the electromagnetic 
component of the plasma (electrons, positrons, and photons) at the moment when 
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neutrino interactions was strong enough so that at the moment of decoupling they 
were in thermally equilibrium state with the same temperature as photons and e . 
According to simple estimates the decoupling temperature, Td ec , for v e is about 2 MeV 
and that for and v T is about 3 MeV. In reality the decoupling is not instantaneous 
and even below T^ ec there is some residual interaction between e ± and neutrinos. 
An important point is that after neutrino decoupling the temperature of the elec- 
tromagnetic component of the plasma became somewhat higher than the neutrino 
temperature. The electromagnetic part of the plasma is heated by the annihilation 
of massive electrons and positrons. This is the well known effect which ultimately 
results in the present day ratio of temperatures, T 7 /T„ = (11/4) 1 / 3 . During primor- 
dial nucleosynthesis the temperature difference between electromagnetic and neutrino 
components of the plasma is small but still non-vanishing. Due to this temperature 
difference the annihilation of the hotter electrons/positrons, e + e~ — *> vv, heats up the 
neutrino component of the plasma and distorts neutrino spectrum. The average neu- 
trino heating under assumption of their equilibrium spectrum was estimated in refs. 
0. [|. Spectrum distortion in Boltzmann approximation was calculated numerically 
in ref. [|J and analytically in ||. In accordance with the latter it takes the form: 



The distortion of the spectra of z/ M and v T is approximately twice weaker. Here 
8f = f-fM. 

An exact numerical treatment of the problem was first done in ref. || and later, 
with a better precision and a corrected expression for the matrix element of one of 
the participating reactions, in ref. The system of coupled kinetic equations (Q) 
governing the neutrino distribution functions with the collision integral of the form 




(5) 



cou 2E 1 ^J 2E 2 (2n) 3 2£ 3 (2tt) 3 2E 4 (2n) 3 
{2^\ Pl + P2 -p 3 - p,)F(h, f 2 , / s , h)S \A\ 2 12 ^ 34 




(6) 
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was solved numerically with the precision about 10 4 . In the expression (|6]) 

F(fu / 2 , /a, U) = fMl ~ /i)(l - / 2 ) - /i/ 2 (l - / 8 )(1 - / 4 ) (7) 

and |^4|i2_>34 is the matrix element squared of the 4-fermion weak interaction. The 
results of ref. [|7| confirmed the shape of spectrum distortion (|5| ). The total relative 
change in neutrino energy densities was found to be 

<Wpo = 0.9%, Sp v ^Jp o = 0A% 1 (8) 

where po is the unperturbed neutrino energy density. 

Naively one would expect that distortion of neutrino energy density at a per 
cent level would result in the similar distortion in the primordial abundances of light 
elements. However this does not occur by the following reason. An excess of neutrinos 
at high energy tail of the spectrum results in excessive destruction of neutrons in the 
reaction: 

n + v e <-> p + e~ (9) 

and an excessive creation in the reaction: 

n + e + <-> p + v (10) 

This nonequilibrium contribution into the second process is more efficient because the 
number density of protons at nucleosynthesis (when T w 0.7 MeV) is 6-7 times larger 
than that of neutrons. So an excess of high energy neutrinos results in an increase 
of the frozen neutron-to proton ratio, r, and in the corresponding increase of A He. 
On the other hand an excess of of neutrinos at low energies results in a decrease of 
r because reaction fliTf ) is suppressed due to threshold effects. It happened that the 
discussed above nonequilibrium spectrum distortion took place in the middle between 
the two extremes and the net influence of these distortion on e.g. 4 He is quite small, 
the change of the mass fraction of 4 He is ~ 10~ 4 . 
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Though quite small, such extra heating of neutrinos may be in principle noticed in 
future high precision measurements of CMB anisotropies || [|. A change in neutrino 
energy density with respect to the standard case would result in a shift of equilibrium 
epoch between matter and radiation, which is imprinted on the form of the angular 
spectrum of fluctuations of CMB. Because of potential observability of distortion of 
neutrino energy density it was recalculated in ref. || where a larger result, than 
found in the previous papers, was obtained. In this connection we |]T0j] repeated our 



calculations with a larger number of integration points, and wider momentum range, 
checked the stability of our calculation procedure and confirmed our previous results 
with the precision of about 10~ 4 . One possible source of disagreement may be an 
incorrect probability of the reactions v a + v a — > v a + v a [a = e, fi, r) used in ref. |§ 
and a smaller number of integration points in the essential region. 

3 Nonequilibrium massive v T and nucleosynthesis 

It is well known that comparison of calculated primordial abundances of light elements 
with observations permits to put a constraint on the expansion/cooling rate at the 
primordial nucleosynthesis (NS) epoch. In particular such arguments allow to limit 
the number of possible neutrino species (or other particles abundant at NS) ]ll], [12], 
13], KM. The present day data seem to exclude one extra neutrino species and possibly 



even 0.3 (for a recent review and analysis see e.g. ref. fllql ). 

Similar arguments permit to put a stringent bound on the mass of v T , considerably 
better than the existing direct experimental limit, vh Vt < 18 MeV |TB|]. Using this 
result and nucleosynthesis data one can conclude that m Vr < 0.5 — 1 MeV, if such 
neutrino is stable at the nucleosynthesis time scale, i.e. t„ t > 100 sec. 

Equilibrium energy density of massive particles is smaller than the energy density 
of massless ones. So if it was the case, then massive v T would effectively correspond 
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to a smaller number of massless neutrinos. However when expansion rate H = a/a 
becomes smaller than the rate of v T - annihilation, equilibrium is no more maintained 
and the number and energy densities of v T become larger than their equilibrium 
values. Because of that one massive v T could correspond to several massless neu- 
trino species. Original calculations of the bound on m UT fTJ§ were made under the 



simplifying assumptions of Boltzmann statistics and kinetic equilibrium of all partic- 
ipating particles. In other words the distribution of massless v e and v T were taken as 
/ = exp(—E/T), while the distribution of massive v T were assumed to have the form: 

f Vr = eM-E/T + (11) 

where the dimensionless (pseudo) chemical potential £ is a function of time only and 
does not depend on the particle momentum. (If lepton asymmetry is vanishingly, 
small the values of £ for particles and antiparticles are the same.) In this approxi- 
mation the complicated integro-differential equations (^H) are reduced to the well 



known ordinary differential equation [18 



h v + 3Hn u = (a ann v)(n^ 2 - nl) (12) 

Here n^ eq ' is the equilibrium number density, v is the velocity of annihilating particles, 
and angular brackets mean thermal averaging. 

The assumption of kinetic equilibrium ( |TT|) is fulfilled if the rate of elastic scat- 
tering at the moment of annihilation freezing, T ann ~ H, is much higher than both 
the expansion rate, H, and the rate of annihilation, T ann . It is generally correct be- 
cause the cross-sections of annihilation and elastic scattering are usually of similar 
magnitudes but the rate of annihilation, T ann ~ a ann n m is suppressed relative to the 
rate of elastic scattering, Y e i ~ <J e ino, due to Boltzmann suppression of the number 
density of massive particles, n m , with respect to that of massless ones, no- However 
in the case of MeV-neutrinos both rates T ann and T e i at the moment of annihilation 
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freezing are of the same order of magnitude. Correspondingly assumption of kinetic 
equilibrium at annihilation freezing is strongly violated. A semi-analytic calculations 
of the deviations from kinetic equilibrium were done in ref. |19|| , where a perturba- 
tive approach was developed. In the case of a momentum-independent amplitude of 
elastic scattering the integro-differential kinetic equation in the Boltzmann limit can 
be reduced to the following differential equation: 

JC" + 2J'C = - A ,V e y ' 2 d v \e~ y d v \eW 2 uyd x fCe^l } (13) 

where x = m/T, y = p/T, prime means differentiation with respect to y, C(x,y) = 
exp(^x' 2 + y' 2 )f m (x, y), and f m is the unknown distribution function of massive par- 
ticles. 

A direct application of perturbation theory (with respect to small deviation from 
equilibrium) to the integro-differential kinetic equation (|2|) is impossible or very diffi- 
cult because the momentum dependence of the anzats for the first order approxima- 
tion to f(p,t) is not known. Numerical solution of exact kinetic equations [|22] shows 
a reasonable agreement with the semi-analytic approach based on eq. fllB"!). 

It can be easily shown that the spectrum of massive v T is softer (colder) than 
the equilibrium one. Indeed if elastic scattering of u T , which would maintain kinetic 
equilibrium is switched-off, the nonrelativistic v T cool down as 1 /a 2 , while relativistic 
particles cool as 1/a, where a(t) is the cosmological scale factor. Since the cross- 
section of annihilation by the weak interactions is proportional to the energy squared 
of the annihilating particles, the annihilation of nonequilibrium v T is less efficient and 
their number density becomes larger than in the equilibrium case. 

After the pioneering paper |17| the frozen number density of v T was calculated 
with increased accuracies in refs. |2(| |21], ^| . The better was the accuracy the larger 
was the calculated frozen number density . In the maximum of n\f> , which occurs 



at m„ T ~ 5 MeV, the difference between the calculations of the papers |17| and is 



almost 50%. So the account of nonequilibrium effects in the distribution of massive 
v T results in a larger frozen number density of v T and in a stronger influence on 
primordial nucleosynthesis. 

Another nonequilibrium effect is an extra cooling of massless v e and due to 
their elastic scattering on colder z/ r , v e tl + v T — > u e ^ + v T . Because of that the inverse 
annihilation v e ^ + v e ^ — > v T + v T is weaker and the frozen number density of v T is 
smaller. But this is a second order effect and is relatively unimportant. 

Considerably more important is an overall heating and modification of the spec- 
trum of v e (and of course of v e ) by the late annihilation v T + v T — > v e + v e (the same 
is true for but electronic neutrinos are more important for nucleosynthesis because 
they directly participate in the reactions (|9| JlO|) governing the frozen n/p-ratio. It is 
analogous to the similar effect originating from e~e + -annihilation, considered in the 
previous section, but significantly more profound. The overfall heating and spectral 
distortion work in the opposite way for m„ T > 1 MeV. An overall increase of the 
number and energy densities of v e and v e result in a smaller temperature of neutron 
freezing and in a decrease of the n/p-ratio. On the other hand a hotter spectrum of v e 
shifts this ratio to a large value, as discussed in the previous section. The latter effect 
was estimated semi-analytically in ref. P3fl , where it was found that e.g. for itl Vt = 20 
MeV the spectral distortion is equivalent to 0.8 extra neutrino flavors for Dirac v T 
and to 0.1 extra neutrino flavors for Majorana v T . The effect of overall heating was 
found to be somewhat more significant p4| , p2| . 

Though the frozen number density of v T obtained in ref. |22| is the largest (in 
comparison to the results of refs. [|Ij], ^0, |2TJ] ) , the influence of nonequilibrium correc- 
tions on nucleosynthesis found in |22j is somewhat weaker than that found in |24], ^1 



in the mass range above 15 MeV. It is possibly related to a larger momentum cut-off 



in numerical calculations of ref. p2| , which gives rise to a smaller neutron freezing 
temperature. For the graphical presentation of the results and comparison with the 
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other papers one can address ref. f22|l . The results of all nonequilibrium calculations 



are systematically and considerably larger than those of the equilibrium ones of ref. 
T7J . These newer and more accurate calculations permit to close the window in the 



mass range 10-20 MeV, which is not excluded by nucleosynthesis if the permitted 
number of extra neutrinos flavors is 1. Now even if 1 extra neutrino is permitted, 
the upper bound on m Vr is about 1 MeV. If 0.3 extra neutrino flavors are allowed, 
the v T mass is bounded from above by 0.3 MeV. These results are valid for the Ma- 
jorana v T . For the Dirac case the mass bound from SN1987 is much more restrictive 
and, moreover the calculations of nucleosynthesis limit on the mass of Dirac u T are 
considerably more complicated because of a larger number of independent unknown 
distribution functions. 

4 Lepton asymmetry and possible spatial variation 
of primordial abundances. 

In the standard nucleosynthesis calculations is usually assumed that neutrinos are not 
degenerate, or in other words, that their chemical potentials are vanishing and their 
distributions are given by the expression ([!]). A justification for this assumption is a 
small value of the baryon asymmetry, but strictly speaking very little is known neither 
from observation nor theoretically about lepton asymmetry. The best observational 
bounds are found from primordial nucleosynthesis (for a recent reference see e.g. ||25|| ). 
Theoretically lepton asymmetry could be as small as the baryon one, especially in 
the models with (B — L)-conservation, but it also may be as large as unity [^, 
p8| . Moreover, the asymmetry could be not only large but also varying by unity at 



astronomically large scales [26 



The recent data [29, 30, pTl B3, B3, am though rather controversial, may possibly 



indicate that the abundance of primordial deuterium changes at the scales of the order 
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of a gygaparsec or several hundred megaparsecs. If the effect is real, there could be 
two possible explanations of it. First, baryon asymmetry of the universe may be not 
a universal constant but a varying function of space points PR I5BI. This possibility 



meets certain problems with the primordial 7 Li-abundance [35| or with the isotropy 
of CMB [[36]. Below we discuss another possible source of a possible variation of 
primordial abundances, namely spatially varying lepton asymmetries |37j. 

It is noteworthy that independently of the data and theory, there is a question 
what is known about light element abundances at large distances. For example what is 
the upper or lower limit on R p , the mass fraction of primordial 4 He, at the distances 
above 100 Mpc? Is R = 50 — 60% or even close to 100% excluded? What is the 
characteristic scale where a large variation of primordial abundances are permitted? 
It is known from observations that the universe is (or better to say was at the early 
stage) very homogeneous energetically. From isotropy of CMB it follows that 

^ < (a few) x 10~ 5 (14) 
P 

A natural implication of the energetical homogeneity is the chemical homogeneity but 
it is not necessarily so. It is interesting to consider a model which gives rise to a small 
cosmological energy variation but to a large chemical variation. 

We assume that chemical potential of neutrinos, especially of u e , are varying on 
the scales above a few hundred Mpc. To explain the possibly observed variation of 
deuterium, the dimensionless chemical potential of electronic neutrinos, £„ e should 
vary by approximately unity. For example £„ e = in our neighborhood and £ Ve — — 1 
in deuterium rich regions. With such variation of electronic asymmetry we immedi- 
ately obtain Sp/p ~ (a few) x 10 -3 , much larger than the bound (||). To save the 
model one has to assume a kind of lepton conspiracy ||37|| , namely if in some space 
region of the universe lepton asymmetry is given by the set of chemical potentials: 

= {<*,P, 7>, (15) 
11 



then in another space region the asymmetry is given by a permutation of a, j3, and 
7. In this case the variation of the cosmological energy density would vanish in the 
first approximation. Though the assumption of lepton conspiracy looks quite 
strong and artificial fine-tuning, it can be rather naturally realized due the flavor 
symmetry, e «-> /i r. 

It can be easily checked that if the variation of 2 H is created by the variation of 
£„ e from to (-1), the corresponding mass fraction of 4 He in deuterium rich regions 
should be larger than 50% []37[ . Such a large variation of helium mass fraction would 
result to a considerable density fluctuations due to different binding energies of helium 
and hydrogen. Rescaling the estimates of ref. |36[ one can find [37J for the fluctuations 
of the CMB temperature : 

ST _ 5 / Rhor 



5 ~nor , . 

T V 10A / v ; 

where A is the wavelength of the fluctuation and Rhor is the present day horizon size. 
The restriction on the amplitude of temperature fluctuations would be satisfied if A > 
200 — 300Mpc//iioo (h 100 = if/100 km/sec/Mpc). Surprisingly direct astrophysical 
effects of such big fluctuations of R p at distances above 100 Mpc cannot be observed 
presently, at least the evident simple ones. 

Another possibly dangerous effect is the differential neutrino heating considered 
in section 2. If chemical potentials of neutrinos are different in different space points, 
their nonequilibrium heating by e + e~-annihilation would also be different. Corre- 
spondingly the photon temperature would also be different. This effect was estimated 
in ref. [j?7|], where it was found that 5T/T k2x 10~ 5 for 5£ v = 1. 

A variation of mass fraction of primordial A He could be observed in the future high 
precision measurements of CMB anisotropics at small angular scales [Q. There are 
two possible effects, first, a slight difference in recombination temperature which log- 
arithmically depends on hydrogen-to-photon ratio, and second, a strong suppression 
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of high multipoles with an increase of R p . The latter is related to the earlier helium 
recombination with respect to hydrogen and correspondingly to a smaller number of 
free electrons at the moment of hydrogen recombination. This in turn results in an 
increase of the mean free path of photons in the primeval plasma and in a stronger 
Silk damping. The position and the magnitude of the first acoustic peak remains 



practically unchanged p8 



This effect seems to be very promising for obtaining a bound on or an observation 
of a possible variation of primordial helium mass fraction. If this is the case then 
the amplitude of high multipoles at different directions on the sky would be quite 
different. The impact of the possible variation of primordial abundances on the 
angular spectrum of CMB anisotropy at low I is more model dependent. It may have 
a peak corresponding to the characteristic scale R > 200 — 300 Mpc or a plateau, 
which would mimic the effect of the hot dark matter. 

5 Neutrino oscillations and nucleosynthesis 

An influence of oscillating neutrinos on nucleosynthesis depends on possible oscilla- 
tion channels. If the oscillations do not create any new neutrino states, and if the 
initial (generated in the early universe) lepton asymmetry is small, the impact of the 
oscillations on nucleosynthesis is negligible. In the case of nonzero lepton asymme- 
try the oscillations between u e , and v T (and their antiparticles) would result in a 
mixing of the different lepton numbers. So that if the oscillations were fast enough 
at nucleosynthesis (NS) and the equilibrium was established, all chemical potentials 
would be equal. If the oscillations at NS were slow, then the asymmetries would not 
be equalized and due to different refraction indices for particles and antiparticles (see 
below) there might be even a significant amplification of asymmetries. 

A more interesting for NS effect takes place if neutrino oscillations produce new 
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neutrino states, e.g. a sterile neutrino (or neutrinos). It may happen if the neutrino 
mass matrix contains both Dirac and Majorana mass terms [|3~9[| . In that case the 
sterile neutrino(s) is (are) just the usual neutrino(s) with a wrong (positive) helicity 
induced by the Dirac mass. If the characteristic time of oscillations is sufficiently 
small, so that thermal equilibrium with respect to formation of new states is fulfilled, 
there would be one or several new neutrino species in the plasma and the only effect 
on NS is the corresponding change in the expansion rate. It is mentioned above that 
one additional neutrino species is forbidden by NS. This condition permits to exclude 



a certain range of the oscillation parameters. In the original treatment of ref. [39| 
the influence of the medium on neutrino oscillations was neglected. In this case the 
characteristic time of oscillations is just the vacuum time: 

E , E/MeV , x 



5m 2 5m 2 /10- 6 eV 2 

The rate of production of new neutrino species is T osc = (r osc sin 2 26^ . If Y osc > H 
then the extra neutrino species would be abundantly produced. 

However for a large and interesting interval of masses and mixing angles the in- 
fluence of the medium cannot be neglected and one should take into account that 



neutrino refraction index in the primeval plasma at NS epoch is not unity |40| : 



n ±-i = M3 xn L^- + c£z- (18) 

E a 

where numerical coefficients Cj are of order unity, Gf is the Fermi coupling con- 
stant, a = 1/137 is the fine structure constant, E is the neutrino energy, T is their 
temperature, and tjl is the leptonic asymmetry of the plasma. There can be differ- 
ent asymmetries for different leptonic charges, then the expression above should be 
correspondingly changed. 

Neutrino oscillations with the account of dispersion effect were considered in refs. 



41| , |42| , [43j , [44j , [45| , (46[. It was shown that the oscillation parameters are roughly 
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speaking bounded by 



sin 4 fl|5m 2 | < 10" 2 eV 2 , if sin 2 # < 0.1 (19) 

and 

\5m\ 2 < 10~ 6 eV 2 , if sin 2 « 1 (20) 



More recent calculations |57, ^8], ^] led to further clarification of the bounds. It was 
shown in particular |48| that spectral distortion of oscillating neutrinos, neglected 
in earlier calculations, is quite essential for accurate determination of the changes in 
primordial abundances due to oscillations. 

A very interesting effect may take place if the MSW-resonance condition is fulfilled 
for oscillations of neutrinos into sterile species. From the expression for the refraction 
index (|18D one can see that the resonance condition is fulfilled either for neutrinos 
or anti-neutrinos depending on the sign of the mass difference. If for example the 
transition of neutrinos into sterile component is enhanced, then the leptonic asymme- 
try in the sector of the usual (not sterile) neutrinos would rise up and the oscillation 
would become more efficient, in turn producing more asymmetry. The equation for 
asymmetry generation has the form 

L = +AL (21) 

where L is the lepton asymmetry and A is a positive coefficient. When the back 
reaction of the oscillation on the initial state of the plasma can be neglected, the 
asymmetry rises up exponentially and can reach the values close to unity. This effect 
was noticed in refs. |5T], and the detailed calculations showing that the effect can 
be quite large was done in ref. ||27|| . If e.g. the asymmetry is generated in electronic 



charge, the impact on primordial nucleosynthesis would be quite significant and in 
particular the limits on oscillation parameters might be less restrictive. Still the 
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mixing angles close to one and relatively large mass differences, 5m 2 > 10 3 , are 
forbidden if less than one extra neutrino species is allowed by NS. 

6 Conclusion 

We see that nonequilibrium neutrino kinetics is quite essential at nucleosynthesis. 
Even for the usual massless neutrinos the deviations from equilibrium are rather 
large, at a per cent level. Though it has very little effect on primordial abundances 
of light elements, about 10~ 4 , the corresponding changes in neutrino energy density 
may be in principle observed in the future high precision measurement of angular 
variation of cosmic microwave background. 

Much more significant are nonequilibrium corrections to the spectra of possibly 
massive tau-neutrinos if their mass lays in MeV region. Exact calculations with all 
nonequilibrium corrections differ from the simpler equilibrium ones as much as by 
50%. The nonequilibrium results are more restrictive and permit to close a window 
for v T mass near 15 MeV which existed in equilibrium calculations if primordial nu- 
cleosynthesis allowed for one extra neutrino species. This permits to put the upper 
bound on Majorana type mass of v T down to approximately 1 MeV. 

Despite the fact that the early universe was very smooth energetically, it is not 
excluded that chemically it is quite inhomogeneous. In particular there is absolutely 
no observational bounds on very large variations of primordial A He at big distances. 
It is definitely worth to obtain from astronomical data any, even very crude limits 
on its mass fraction, R p . Again CMB measurements could be very helpful in this 
respect. Varying R p would give rise to a different amplitudes of high multipoles at 
different directions on the sky. 

In the case of oscillations of the known neutrinos between themselves primordial 
nucleosynthesis does not permit to put any interesting bound on the parameters of the 
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oscillations. However if neutrino may oscillate into new (sterile) ones, large mixing 
angles and large mass differences are excluded. 
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